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Abstract. The brain is one of the most studied and highly complex sys- 
tems in the biological world. It is the information center behind all ver- 
tebrate and most invertebrate life, and thus has become a major focus 
in current research. While many of these studies have concentrated on 
studying the brain directly, our focus is the structure of the brain itself: 
at its core an interconnected network of nodes (neurons). A better under- 
standing of the structural aspects of the brain should elucidate some of 
its functional properties. In this paper we analyze the brain of the nema- 
tode Caenorhabditis elegans. Consisting of only 302 neurons, it is one 
of the better-understood neural networks. Using a Laplacian matrix of 
the 279-neuron "giant component" of the network, we use an eigenvalue 
counting function to look for fractal-like self similarity. This matrix rep- 
resentation is also used to plot (in eigenfunction coordinates) both 2 and 
3 dimensional visualizations of the neural network. Further analysis ex- 
amines the small-world properties of the system, including average path 
length and clustering coefficient. We then test for localization of eigen- 
functions, using graph energy and spacial variance. To better understand 
these results, all of these calculations are also performed on random net- 
works, branching trees, and known fractals, as well as fractals which have 
been "rewired" to have small-world properties. This analysis is one of 
many stepping-stones in the research of neural networks. While many of 
the structures and functions within the brain are known, understanding 
how the two interact is also important. A firmer grasp on the structural 
properties of the neural network is a key step in this process. 
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Author Summary 

The brain is the biological center driving all animal life, and certainly 
our own. Though fundamentally an interconnected system of nodes, the 
information-transfer that occurs inside of these networks perpetuates life it- 
self. As a result, the brain is the focus of immeasurable scientific research. 
While many studies clinically analyze the brain and its function, our goal 
is to look at the brain on a simpler level- as a network of connected points. 
Studying the structural connectivity and network properties within the brain 
will help further our understanding of the organization underlying its many 
functions. In this paper, we study the neural network of the nematode worm 
Caenorhabditis elegans. It is a well-understood system consisting of 302 
neurons, making it an excellent candidate for our research. (The human 
brain is composed of billions of neurons, the connections between which 
are not completely known). Through the course of our study, we look for 
both self- similarity and small-world characteristics in the structure of the 
C. elegans neural network. A better understanding of the brain's physical 
configuration will help to reveal the mechanisms behind its structure and 
function. 



1. Introduction 

Fractal theory has become an increasingly popular topic of both debate 
and research in recent years. Beginning with Mandelbrot's discussion of 
Britain's immeasurable coastline [1], fractal theory has found applications 
in both the mathematics and scientific communities. In the geometric sense, 
fractals are objects that exhibit self-symmetry: they exhibit the same pattern 
on increasingly smaller scales. In other words, magnifications of smaller 
portions resemble the whole object. 

More recently, fractal theory has found applications in the biological realm. 
Kinetics of ion channels have been modeled with fractal structures [|2l[3]|. 
Fractal dimension has been used to analyze human EEG signals [|4| as well 
as the complex morphology of living cells [|5l|6||. The applications of fractal 
theory in neuroscience have been a particularly prevalent topic of research 
[HI m |9l [ini [m . OUal cells have been analyzed in-depth using fractal di- 
mensions and modeling [rT2l[T3l[T4ll . Dendritic branching has been shown 
to have self- similarity [[T5l[T6]| . and other studies have examined fractal pat- 
terns in neuron connectivity fVl]. Further applied research has used three- 
dimensional fractal structures to approximate the white matter surface of the 
human brain, based on MRI images IITSll . 
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In this paper we use a graph-theoretical approach to probe the structure of 
the Caenorhabditis elegans neural network for self-similar structures. Ad- 
vances in graph theory have proven highly useful in analyzing complex neu- 
ral networks [Il9l|20ll2D as well as underlying motifs in the brain [I22ll23ll . 
These methods have made considerable contributions to our understanding 
of the structure and function of one of biology's most intricate and important 
systems. In this paper we apply these techniques to a physical map of the 
C. elegans brain. With a well-connected component of only 279 neurons, it 
is an excellent candidate for graph theoretical research on a complete-brain 
model. While [24J presents a geometric structure of this nematode brain, our 
research continues that of [25J in which Chklovskii et al. propose a finalized 
schematic of the C. elegans neural network. 

The C. elegans brain is composed of three types of neurons: sensory neu- 
rons, motor neurons, and interneurons. Two types of connections exist be- 
tween these neurons: chemical synapses and gap junctions. The gap junction 
network, which sends electrical signals via ion transport, is an undirected 
system. Conversely, chemical synapses possess clear directionality [25] . 
We are only interested in studying the overall connectivity between neurons. 
Thus in order to analyze the structure of the C. elegans neural network, we 
consider only the skeleton of the brain's organization. Although some neu- 
rons share multiple points of contact (they have a multiplicity of connection 
> 1) we consider this a single connection. While chemical synapses send 
directional signals, we only observe that two neurons are connected: re- 
gardless of that connection's direction. As a result, we are able to study a 
weakly connected network representing only the framework of connections. 
(See Methods) 

In order to index each of these connections we used the graph Laplacian 
matrix, L = [kj]. For a graph, G, we define as the degree of a vertex v: 
the number of total connections. (Note that each vertex represents a neuron). 
If vertex u is connected to vertex w then lu,w = — 1 and = — 1. These 
correspond to the entries in the u^^ row and w^^ column, and the w*^ row 
and u^^ column. Furthermore, L^^^ = d^. All other entries of matrix L 
are 0. We also generated similar matrices for randomly generated networks, 
random branching trees, and known fractals. 

The original goal of this study was to examine the structure of the C. ele- 
gans neural network for self-similarity. Through the course of our research 
we analyzed many other network properties, and compared these results with 
similar calculations on other systems. Although fractal theory has repeat- 
edly been applied in neuroscience, in studying the structure of the nematode 
brain the results are not as simple as saying "fractal" or "not- fractal." Instead 
we search for self-similar structures, as there is no precise definition of what 
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it means for a network to have "fractal properties." Nevertheless, we uncover 
several interesting properties of the C. elegans neural network. 

2. Results and Discussion 

Before we could begin our analysis, we had to construct a variety of Lapla- 
cian matrices. For our C. elegans model, we derived a Laplacian matrix 
from the adjacency matrices used in [25J (See Methods). We already had 
the Laplacian matrices of known fractals such as the Sierpinski Carpet, the 
Octagasket, the Hexacarpet, and the Sierpinski Gasket (specific subdivisions 
were chosen which had similar numbers of vertices as the C. elegans neural 
network). Next we wrote a series of MATLAB programs which would gen- 
erate Laplacian matrices with specific conditions. One produced that of a 
random network when given a defined number of vertices (n) and probabil- 
ity of connection (p). Another produced the Laplacian matrix of a random- 
branching tree, given the total number of vertices (n) as well as the maxi- 
mum number of branches possible from any point (m). One final program 
randomly rewired these networks, moving connections from one vertex to 
another with a given probability, p. Details on these programs can be found 
in Methods O and El 

2.1. The Eigenvalue Counting Function. To begin our analysis, we ap- 
plied the eigenvalue counting function to our Laplacian matrices. The eigen- 
value counting function is a cumulative distribution function on the spec- 



trum of a matrix (see Methods 3.2 1. It computes the set of all eigenvalues of 
a given matrix (the spectrum), and counts the total number of eigenvalues 
less than the given input. Plotting this function gives an expedient way to an- 
alyze how the graph of a given matrix should be generally organized, based 
on the spectrum of the graph Laplacian ^26\. Figure |2]T shows the plots of 



the eigenvalue counting function on a variety of our Laplacian matrices. 

Upon examination of these graphs, a few patterns immediately emerge. 
The first is the presence of step-like portions of those graphs corresponding 
to known fractals (Fig. l(b)[|l(d)]|l(f)[[T(h)| ). These sections of slope-zero 



correspond to gaps in the spectrum of the graph. This feature is closely 
linked to self-similar fractal structures [|27| . Although the eigenvalue count- 



ing function plot of the C. elegans brain (Fig. 1(a) ) does not show these de- 
finitive spectral gaps, this does not conclusively eliminate the possibility of 
self-similar structures existing within the neural network. At the same time, 
however, this does indicate that the nematode brain is not strictly fractal- like 
in structure. 

Another interesting pattern appears in the graph corresponding to the 



random-branching tree (Fig 1(c)). There is a large vertical jump at x = 1, 
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(a) C. elegans neural network 




(b) Sierpinski Gasket, Level 5 










(c) Random Branching 

n = 279, m = 10 


Tree 


(d) Hexacarpet Level 3 










(e) Random Network n = 279, p = 0.07 


(f) Octagasket, Level 3 



















(g) Sierpinski Gasket Rewiring p — 0.15 (h) Sierpinski Carpet, Level 3 



Figure 2.1. The Eigenvalue Counting Function 



with a change on the y-axis of approximately 200. This indicates that for this 
graph Laplacian, the eigenvalue 1 occurs with extremely high multiplicity. 
This is caused by the nature of the tree's structure. With a finite number of 
points (in this case 279) a large number of these points are endpoints: ver- 
tices at which no further branching occurs. These points are only connected 
to one other: their "parent" vertex. The lack of this spectral "jump" in the C. 
elegans brain indicates that there are very few singly-connected endpoints 
in the system. Rather, the neural network is highly inter-connected, with 
structural properties baring little resemblance to the branches of a tree. 
On the other hand, the eigenvalue counting patterns of the C. elegans 



neural network do resemble those of the random network (Fig. 1(e)) and 



the rewired Sierpinski Gasket (Fig 1(g) I. While a similarity in these patterns 
(as opposed to the aforementioned dissimilarities) cannot conclusively point 
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to similar structural properties, these eigenvalue counting functions are a 
prerequisite to the next step in our analysis, Weyl ratios. 



2.2. Weyl Ratios. In order to find Weyl ratios for our network, we plotted 
the respective eigenvalue counting functions on a log-log scale. We then 
searched for an interval on which this graph is roughly linear. We determined 
a line of best fit for this section, and found its slope, a. We then used this a 
to produce a Weyl ratio for each graph (see Methods 3.2[ ). A Weyl ratio is 



essentially a rescaling of points in our eigenvalue counting graph, providing 



a more revealing visualization. Figure 2.2 shows the Weyl ratios for each 
of our Laplacian matrices. Periodicity in the Weyl ratio plot suggests the 
presence of self-similar structures. For more on Weyl ratio analysis of known 
fractals, see fi28l . 



As expected, the Weyl ratios of known self-similar fractals (Fig. 2(b)[ 



|2(d)[ |2(f)[ 2(h)[ ) show a high degree of organization. That of the Sierpinski 



gasket in particular (Fig. 2(b)), shows unmistakable periodicity. Although 
not as readily noticeable, clear patterns exist in the Weyl ratios of other frac- 
tals as well. 

It should be noted that the Weyl ratio graph of the branching-tree (Fig. 



2(c)) is difl'erent from those of other networks. Although the tree structure 
lacks the high ordering of known fractals, there is a clear distinction between 
this special case and all other networks: what we will call "looping". In 
fractals, random networks, and the C. elegans brain alike, the vertices are 
highly interconnected. Many cyclic paths exist which allow a signal to arrive 
back at a starting vertex by traveling through a series of other vertices. Trees, 
on the other hand, lack this feature: only one path exists between any two 
points. Not only does this create a unique Weyl ratio pattern, but it also 
suggests that neural networks are not strictly branching structures, as one 
might expect. 

While the Weyl ratio pattern of the C. elegans brain ( Fig. |2(a)| ) is dissim- 
ilar from that of a branching tree, it also diff"ers greatly from those of known 
fractals. While several cases of slight periodicity could be argued for, this 
evidence is not definitive enough to indicate concrete self-similarity in the 
C. elegans neural network. However, there does exist similarity between the 
Weyl ratio patterns of the C. elegans neural network, the randomly generated 



network (Fig. 2(e) ), and the rewiring of the Sierpinski Gasket (Figure 2(g) ). 
(It should be noted that many random networks and SG rewiring s were gen- 
erated, these two examples were chosen as representatives). Although the 
significance of examining a "rewired" fractal structure will be explained in 



Section 2.4 of this paper, it is important to note that this likeness in Weyl 
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(a) C. elegans neural network (b) Sierpinski Gasket, Level 5 




(c) Random-Branching Tree (d) Hexacarpet Level 3 

n = 279, m = 10 




(e) Random Network n = 279, p = 0.07 (f) Octagasket, Level 3 




(g) Sierpinski Gasket Rewiring p = 0.15 (h) Sierpinski Carpet, Level 3 

Figure 2.2. Weyl Ratios 



ratio patterns can suggest a structural similarity. While it seems counterin- 
tuitive that the nematode brain would be a random arrangement of neurons, 
the evidence found in the Weyl ratios alone cannot conclude that this is (or 
is not) indeed the case. Our further analysis will make distinctions between 
these 3 networks. 
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2.3. The Eigen-Projection Method. Next, we wanted a way to visualize 
the neural network. We did this by embedding a graph of all neurons in 
Euclidean space via the eigen-projection method explained in [l29ll . similar to 
those processes described in l^30ll3T1l . This spectral approach to visualizing 
graphs utilizes the eigenfunctions of degree-normalized Laplacian matrices 



(see Methods 3.3). 



The eigen-projection method plots the vertices of a graph using the eigen- 
functions of its Laplacian matrix as a coordinate system (it is also referred 
to as 'plotting in eigenvector coordinates'). This method essentially projects 
the graphs into a smaller Euclidean space, typically either or M^, using 



appropriate eigenfunctions. See Methods 3.4 for a more rigorous descrip- 
tion. 

After embedding each vertex in either 2- or 3-dimensional space, we then 
represented neuronal connections (or network connections) with line seg- 
ments between the appropriate points. In the case of our C. elegans brain 
diagram, we also used the same color-coding as Il25l : where red represents 
sensory neurons, green are motor neurons, and blue indicates interneurons. 
Lastly, we labeled the points with the corresponding neuron name abbrevi- 
ations. This was done using a slight variation of the VISUALIZE program 
used by Chklovskii and Varshney, available at [[32l|. 



The eigen-projection visualizations (Figure 2.3 ) allow us to make further 



distinctions between the C. elegans brain and other networks. In support of 



our previous observations, it is again clear that the nematode brain (Fig. 3(a) 



and 3(b) is not strictly fractal in structure. The eigenfunction graphs of the 



Sierpinski gasket (Fig. |3(c) and 3(d) ), once again display characteristics ex- 



pected of self-similar fractals: a high degree of ordering and self- symmetry. 

While, the eigenvalue counting function and Weyl ratios showed little dis- 
tinction between the C. elegans brain, and a random graph, these eigen- 
projections make a clear differentiation between the two. The random graph 
(Fig. 3(e) and |3(f)[ ) appears, as expected, more or less a scatter of points. The 



C. elegans brain, however, shows a definite structure with organized connec- 
tivity. Thus while our previous results appeared more or less inconclusive, 
these eigen-projection techniques suggest that the structure of the C. elegans 
neural network is not a randomly connected system of neurons. On the other 
hand, the C. elegans neural network maintains its resemblance to a rewired 
Sierpinski gasket when plotted in eigenfunction coordinates (Fig. |3(g) and 



3(h)[ ). While there is no effective way to quantify this heuristic similarity, it 



sustains its interest experimentally. This additional similarity continues to 
suggest the presence of structural parallels. 

The eigen-projection method allows us to view not only the structural or- 
dering, but also the functional organization of the C. elegans neural network. 
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(a) C. efegani neural network, ((^2, '/'s) (b) C. elegans neural network, 



(c) Sierpinski Gasket, Level 5, ((^2, (y^s) (d) Sierpinski Gasket, Level 5, 



(e) Random Network n = 279, p = 0.07, (f) Random Network n = 279, p = 0.07, 

('^52,<y33) (¥'2,'^53,</54) 



(g) Sierpinski Gasket Rewiring p = 0.15, (h) Sierpinski Gasket Rewiring p = 0.15, 



Figure 2.3. The Eigen-Projection Method 



It is clear that the neurons are arranged roughly by neuron type. There is a 
distinctive cluster of motor neurons (green), a larger sub-component of sen- 
sory neurons (red), and interneurons interspersed throughout the network 
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(blue). This indicates that although the brain may not posses the strict self- 
similarity of a fractal structure, it is indeed highly organized. The neuron 
grouping evident in our visualization has most likely developed for entirely 
functional purposes: producing more efficient signal transfer and more or- 
ganized communication between neurons. 



2.4. Small-World Network Properties. When analyzing graphs as net- 
works of nodes, it is useful to consider the nature of the graph's connections. 
To do this, we consider two functions defined on graphs: average clustering 
coefficient and average path length. The clustering coefficient, Cy, of a vertex 
V, is the probability that any two vertices neighboring v are also connected 
to each other. Formally, this can be calculated by counting the number of 
edges between neighbors of v and dividing by the number of total possible 
connections that could exist between neighbors of v. (For more details, see 
Methods |3.5| ). For a graph G, the average clustering coefficient, c, is the 
average of Cy over all vertices. 

The path length between two vertices u and v is the shortest path along the 
graph's edges connecting u and v (Note that this path usually travels through 
a number of other vertices). Using Djisktra's algorithm, it is possible to 
rigorously determine the shortest path between a given vertex and each other 
vertex on the graph. By repeating the algorithm for each node on the graph, 
it is possible to determine the shortest path between each pair of vertices. It 
follows naturally that the average path length, /, is calculated by finding the 
arithmetic mean of the shortest paths between each pair of vertices on the 
graph. We use both c and / to analyze small-world behavior. 

One prominent theme in modern graph-theoretic and fractal research is 
the "small- world phenomenon". This phenomenon is best described as the 
tendency for certain networks to have a much shorter path length than intu- 
ition suggests. For example, the term "Six Degrees of Separation" suggests 
that any two people on earth are no more than six "steps" away from each 
other, even though it seems that this number should be much larger than six. 
This phenomenon can be explained by the generalization that humans live 
in relatively large and tight-knit social networks, and thus have a relatively 
high c value. However, each person also likely possesses a few long dis- 
tance relationships, allowing / to attain very small values without lowering 
the value of c significantly. Thus, small-world networks are (generally) de- 
fined as networks which have a much higher c value than random networks, 
but maintain a value of / only slightly larger than that of a random network 



Small-world networks arise quite often in the natural sciences, as they al- 
low for the efficient transfer of information while maintaining a certain level 
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Graph 
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Average Path Length 


Sierpinski Gasket, Level 5 


0.4495 


17.3721 


Random(Sierpinski Gasket) 


0.0104 


5.748 


Sierpinski Gasket Rewire p = 0.15 


0.2843 


7.3833 


Random(SG Rewire) 


0.0104 


5.748 


C. elegans Neural Network 


0.3371 


2.5377 


Random(C. elegans Neural Network) 


0.0581 


2.3458 



Figure 2.4. Clustering Coefficient and Path Length 



of complexity. There is a great deal of research which suggests that neu- 
ral networks possess small- world properties [33l|34l. Our findings propose 
that C. elegans is no exception. As Figure |24] shows, the C. elegans neural 
network has an average path length only slightly larger than that of its asso- 
ciated random network (see Methods |3.6| for how these 'associated random 
networks' were constructed). At the same time, the average clustering co- 
efficient for C. elegans is six times larger than that of its associated random 
network. This being the case, the neural network of C. elegans satisfies the 
small- world properties as defined by [35 j. 

It should be noted that this discussion of small-world networks is relevant 
to our discussion of self-similar fractals. Certain research [36J suggests that 
there is an apparent dichotomy between fractal structures and small-world 



networks. To illustrate this idea, refer again to Figure 2.4 It is apparent 
that the S*'* level Sierpinski Gasket is not small-world in nature: its average 
path length is significantly larger than that of its associated random network. 
Before we continue our discussion of small- world networks, we must define 
the neighborhood of a graph. 

The neighborhood of a graph, H{m) , where m is a positive integer, is also 
useful in analyzing small- world networks. For our purposes a neighborhood 
of size m around a vertex v is the set of all vertices that can be reached in m 
steps or less from v (see Methods [3 .7 [ for a rigorous definition). The behavior 
of H{m) can tell us quite a bit about structure of a graph, G. If we look at 
the growth characteristics of H{m) for small m, we begin to understand the 
localized neighborhoods of G. Certain research suggests that H{m) grows 
exponentially for small m if G is a small- world network [[36l. Likewise, 
H{m) tends to grow polynomially if G is a self-similar fractal. Figure [5(a) 



shows a plot of H{m) for the C. elegans neural network. Figures 5(b)| and 
5(c)| show the same plot on log-log axes and linear- log axes, respectively. 



Unfortunately, the number of relevant points is simply too small to make 
a concrete statement as to whether or not this H{m) shows exponential or 
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5 6 7 



(a) C. elegans neural network (b) C. elegans neural network, log-log 



5 10 15 20 35 30 35 



(c) C. elegans neural network, linear-log (d) Sierpinski Gasket, Level 5 

Figure 2.5. Graph Neighborhoods 
polynomial growth. On the other hand, H{m) is clearly polynomial for the 



Sierpinski Gasket (Fig. 5(d)). In fact, it appears roughly linear. 



Thus far, our research suggests that C. elegans neural network does not ap- 
pear to structurally resemble a tree, a random graph, or a self-similar fractal. 
However, it does appear to possess small- world properties. This realization 
motivated our work with network-rewiring, related to that done by Watts and 
Strogatz. In [35J they showed that moving connections in an ordered net- 
work, with a certain probability p, led to some interesting changes in graph 
structure. Namely, when p is small, a slight increase in p causes a large drop 
in / but does not change c appreciably: thus the network takes on small- world 
characteristics. Intuitively this can be explained by the fact that these sparse 
random connections don't change a graph's strong localized structure, but it 
is now easier to travel long distances via these new connections which can 
span large gaps. 

This was our motivation in analyzing the rewired Sierpinski Gasket (seen 
throughout previous sections). Self-similar fractals are highly ordered net- 
works, as were the ordered graphs analyzed in OSll . By rewiring certain 
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known fractal structures, we were able to observe if the resulting small- 
world network was comparable to our C. elegans neural network, which has 



been demonstrated to have small-world properties (See Methods 3.8 for our 
rewiring algorithm). This would suggest that the C. elegans brain has lo- 
calized self- similarity with interspersed gap-spanning connections, allowing 



for more efficient signal transfer. Figure 2.4 shows the results of rewiring the 
^th Sierpinski Gasket. Note that the random graphs associated with the 
Sierpinski Gasket and its rewiring are identical because rewiring preserves 
the number of connections, and thus the average degree. It is clear that the 
rewired Sierpinski Gasket is an example of a small- world network based on 
our definitions, and thus maintains interest as a comparison to the C. elegans 
neural network. While previous evaluations inferred similarities between 
these two networks, we now know that both show small- world properties as 
well. 

2.5. Energies and Spacial Variances. Next, we analyzed energies and spa- 
cial variances on each of our graphs. Using an eigenfunction of a graph's 
Laplacian, Lp, one can calculate a graph energy specific to Lp (See Methods 



3.9 ). Knowing the resistance between any two vertices and a constant 7, one 



can also use this (y9 to calculate a spacial variance (See Methods 3.10). These 
two quantities allow us to observe localization of eigenfunctions. This lo- 
calization occurs when eigenfunctions are approximately zero except for in 
a localized region. Localized eigenfunctions are a feature of certain self- 
similar fractals, such as the Sierpinski Gasket [37J. 

The energies and spacial variances (at 7 = 1) of the eigenfunctions of the 
unnormalized Laplacian of the C. elegans neural network were compared to 
those of other graphs. These included random graphs, fractal graphs, small- 
world networks, and random trees. Distributions of these values were plotted 
for each network. Any resemblance in the distributions of the energies and 
spacial variances could indicate some similarity in the structure of the graphs 

m. 



The energies and spacial variances of the random graph (Fig. 6(c) and 
|6(d)[ ) were quite diff"erent from those of the neural network. In the random 
network, the distributions of these values tend to be peaked and symmetric, 
whereas those of the neural network are more spread-out and skewed (Fig. 



6(a) and 6(b) ). This evidence suggests that the neural network of C. elegans 
is not randomly distributed. Intuitively, this is what one would expect, and 
is consistent with our earlier observations. 

Next, we continued our comparison of the neural network of C. elegans 
to known self-similar fractals, once again using the 5^^ subdivision of the 
Sierpinski Gasket. The energies and spacial variances of this fractal (Fig. 
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ft. 




(a) C. elegans neural network energies (b) C. elegans neural network spacial vari- 
ances) 





(c) Random graph (n = 279, p = 0.07) (d) Random graph (n = 279, p = 0.07) 
energies spacial variances 




JJ 



(e) Sierpinski Gasket, Level 5 energies (f) Sierpinski Gasket, Level 5 spacial vari- 
ances 



I 




(g) Random Tree {n = 279, m = 10) en- (h) Random Tree (n = 279, m = 10) spa- 
ergies cial variances 





(i) Sierpinski Gasket Rewiring (p = (j) Sierpinski Gasket Rewiring (p 
0.15), energies 0.15), spacial variances 



Figure 2.6. Graph Energies and Spacial Variances 



6(e) and 6(f) ) vary considerably from those of the neural network. Although 



the spacial variances of all graphs were roughly symmetric, the variances of 
the neural network were much lower than those of the Sierpinski Gasket. 
On the neural network, the spacial variances were all less that 1 x 10^'^, 
whereas those of the fractal network ranged between 1 x 10^^ and 8 x 10^'^. 
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Furthermore, the distribution of energies of the neural network was skewed 
right and ranged from to almost 100, while the energies of the Sierpinski 
Gasket were very jagged (no resemblance to a normal distribution) and con- 
tained many gaps. In addition, all of these energies were in a narrow range 
between and 6. These differences are in agreement with earlier results, 
suggesting that the nematode brain is not strictly fractal in structure. 

We also compared the C. elegans neural network to a random branching 
tree (Fig 6(g) and |6(h)[ ). The energies and spacial variances of the random 
tree are again different from those of the neural network. The energy distri- 
bution of the tree was peaked, with most values falling between and 10, 
whereas those of the neural network were distributed over a wider range, 
from to 100. The energies of the neural network are skewed right. The 
differences in spacial variances were even greater. The spacial variances of 
the neural network lie between and 1 x 10^^, whereas those of the random 
trees are spread over a larger range, up to about 1.5 x 10^^. The distribution 
of the neural network is again skewed right, whereas that of the random tree 
is very much skewed left. Once again, the neural network does not exhibit 
the same structural properties as a random-branching tree. 

Thus far, the "small-world" rewired Sierpinski gasket has shown the most 
similarities to the C. elegans neural network. Once again we analyzed the 
rewiring ofthe fifth level Sierpinski triangle, with p = 0.15. This probability 
was chosen to generate a small-world network without disturbing much of 
the self-similar structure. After rewiring, the energies of the graph (Fig. 



6(i) I were distributed in nearly the same pattern as those of the original self- 



similar fractal, only the values were shifted slightly to the right. The same 
is true of the spacial variances (Fig. |6(j)[ ), except that this distribution was 
shifted to the left. Although the range of the distributions on the rewired 
fractal is closer to that of the neural network, the distributions themselves 
remain markedly different from those of the C. elegans brain. This leads 
us to believe that the structure of the C. elegans neural network, although 
small-world in nature, does not structurally resemble the Sierpinski Gasket. 
On the other hand, the spacial variances of the C. elegans neural network 
are considerably lower in magnitude than those of the Sierpinski Gasket. 
This suggests localization of eigenfunctions, which could in turn indicate 
the presence of self- similarity in the network. 



2.6. Conclusions. Using a variety of mathematical techniques including 
the eigenvalue counting function, Weyl ratios, and the eigen-projection method, 
analyzing small- world properties, and calculating of graph energies and spa- 
cial variances, we were able to uncover some structural properties of the C. 
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elegans neural network. Although much previous research has been dedi- 
cated to applying fractal theory to neuroscience, our results suggest other- 
wise. The known structure of the C. elegans neural network does not exhibit 
properties characteristic of strictly self-similar fractal networks. While some 
evidence suggests that there may be localized instances of self- similarity, we 
cannot quantify or definitively state such a conclusion. Our research found 
additional structures which the C. elegans neural network does not resemble: 
it does not exhibit the branching properties of a tree, nor does it constitute a 
randomly connected network. 

As to what we can conclude, our research is consistent with related work, 
showing that the brain exhibits small-world network properties. Further- 
more, the network has highly localized eigenf unctions, which could suggest 
the presence of self-similar structures. Further research would be required 
to determine the nature of these localized eigenfunctions. Although the C. 
elegans neural network does not appear to be random, tree-like, or fractal in 
structure, it is certainly highly ordered. This organization most likely aids 
in functional efficiency of the system. Further research is needed to deter- 
mine a more refined view of the brain's structural properties. Although C. 
elegans has proven to be a useful model organism, with a well-defined map 
of its neural network, this network consists of only 279 nodes. While this 
makes the system fairly efficient to study computationally, this small number 
of nodes makes network-analysis somewhat limited and rather unrevealing. 
However, due to the difficulty in determining the exact layout of each neuron 
in the network, very few consistent complete-brain maps exist at this time. 
Similar analysis as that provided in this paper, applied to a more complex 
or higher-order neural network, could potentially show more conclusive re- 
sults. 



3. Methods 

In order to analyze only the framework of the C. elegans neural network, 
we constructed a Laplacian matrix derived form the adjacency matrices in 
Il32ll . We wanted to look only at connections between neurons, regardless 
of type or direction. The network of chemical synapses sends signals in one 
direction only, resulting in a non-symmetric adjacency matrix, C. In order 
to disregard this directionality, we added this matrix to its own transpose, C, 
creating a symmetric matrix indexing all chemical connections. We added 
this matrix to the adjacency matrix of the gap junction system, G (which is 
already symmetric because these connections are bidirectional). 



B = [bij] = {C + C') + G such that i,j < 279 
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We then normalized all non-zero entries of this combined matrix, B, to 
be 1 (in order to avoid multiplicity of connection), resulting in an adjacency 
matrix representing only the framework of the entire network, A. 

A — [uij] where aij — hj/bij — lif bij > 0, otherwise a^j- — 0( when bij 
It is then simple to produce a Laplacian matrix, L, as shown below: 

j=279 

dj = ^ ttij for each j < 279 

1=1 

Note that dj is the degree of each vertex j. The degree matrix, D, is now 
defined as: 

D — [dij] where dij — dj when i — j otherwise dij — 0(when i ^ j) 
Then the Laplaican matrix, L, is given by: 

D-A 

3.1. Generating Random Graphs and Trees. Throughout this project, we 
compared the C. elegans neural network to both random graphs and branch- 
ing trees. In order to generate a Laplacian matrix representation of the ran- 
dom graphs, we used the following algorithm: 

First, we fixed the number of vertices, n, and the probability of connec- 
tion, p. We then constructed an empty n x n matrix, R = [rij]. 

For each r^j such that i < j we assign a random value aij such that 

< Oij < 1 for all i,j<n.lf ttij < p then rjj = 1, otherwise r^j = 0. 
To produce an adjacency matrix of this graph. A, we must add R to its 

own transpose: 

A = R + R' 

Using this adjacency matrix we can construct a Laplacian matrix using the 
method described previously. 

The algorithm for producing the Laplacian matrix of a random-branching 
tree is more involved. Again, we fix the number of vertices, n, and also 
specify the maximum number of "children" from any given branch-point, 
m. We create an empty n x n matrix, T — [tij] 

We begin by generating a random integer ai such that < ai < m, and 
take ti^i = tti. This corresponds to the first vertex having |ai| branches. To 
represent these branches in the matrix, we take iij = — Ifor j = 2, • • • ,ai+ 

1 and ti^i — —1 for i — 2, ■ ■ ■ , oi -|- 1. 
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Next we move to all subsequent vertices. Because no "looping" exists in 
the structure of the tree, each node can only be connected to its parent vertex 
and its "children" vertices. We take S = {j : Uj = for all i < n} Then k, 
where k =min(S') is the smallest-labeled node which does not have a parent 
vertex, i.e. the first column with all entries corresponds to the first point 
not yet connected. (Note in the case of vertex 2, /c = oi + 2). This vertex k 
is the first "offspring" from the next branch-point. 

Now, as above, for each remaining vertex v we choose another random 
integer, a^, such that < < min(m, n — k + 1) and take t^^y = + 1. 
(Note that vertex v has children, however a^, + 1 is the degree of node 
V, taking into account its parent-connection). To represent the "offspring" 
branches of this vertex v, we use the following formula: 

ti^y = —1 for i = k, k + 1, ■ ■ ■ ,k + {ay — l) 

and 

tyj = —1 for j = k, k + 1, ■ ■ ■ ,k + (ay~l) 

We use min(m, n — k + 1) when choosing to avoid adding more vertices 
than the n which we originally fixed. This algorithm, when repeated for each 
vertex v, produces the Laplacian matrix of a random branching tree. 

3.2. The Eigenvalue Counting Function and Weyl Ratios. For a given 
graph Laplacian matrix, L, the eigenvalue counting function, N{x) is a cu- 
mulative frequency function on the spectrum of the matrix where: 

N{x) = #{Aj < x} where each Xj is an eigenvalue of L 

The growth of N{x) is approximately x", thus the relevant portion of each 
graph, when using a logarithmic scale, appears to be linear. A line of best 
fit was found for each relevant interval, and the slope, a, calculated. Using 
this a, we plotted the Weyl ratio, W{x), such that: 

W{x) = N{x)/x'^. 

These Weyl ratios allow us to examine the spectrum of each matrix, looking 
for elements such as symmetry and periodicity. ||28]| 

3.3. Normalizing a Laplacian Matrix. We used two different forms of 
the graph-Laplacian matrix: the standard Laplacian matrix and the degree- 
normalized Laplacian. In the case of eigen-projections, we utilize the degree- 
normalized matrix. We define the degree matrix, D, as before: a diagonal 
matrix whose non-diagonal elements are 0, and each entry djj is the degree 
of the j^^ vertex. Using both the standard Laplacian, L and its corresponding 
degree matrix, D, we produce the degree-normalized Laplacian, Q: 

Q = D-^'^LD-^/^ 
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(a) Un-normalized Laplacian, ((/?2,¥'3) (b) Normalized Laplacian {lp2t^z) 



Figure 3.1. C. e/egani' neural network 



Using the degree-normalized Laplacian has many aesthetic advantages, 
as shown in Figure [T(a)] and [T(b)| The normalized matrix also has all eigen- 
values A. such that < A, < 2. 



3.4. Graphing in Eigenfunction Coordinates. We found all eigenvalues, 
Afc, and their corresponding eigenfunctions, tp^, for each matrix. Given two 
eigenfunctions and ^pj, (such that i ^ j) we then plotted the ordered 
pair {(pi{n), Pj{n)) for each n from 1 to 279, as described in [|38ll . The first 
eigenvalue of any Laplacian matrix is always 0, corresponding to a constant 
eigenfunction. Thus we only consider pi and pj with i,j > 2. Edges were 
then added between points to represent relevant connections, and the same 
color-coding as f25\ was used: where red represents sensory neurons, green 
are motor neurons, and blue indicates interneurons. The same process was 
then repeated in three dimensions, plotting ((pi(n), Pj{n), pk{n)) for some 
i, j, k >2, such that i j ^ k. 



3.5. Clustering Coefficient. The clustering coefficient is a common mea- 
sure for vertices on a graph. It is typically measured on graphs with un- 
weighted edges. For a graph G and a given vertex v, let denote the number 
of connections that exist between the neighbors of v. Take as the number 
of neighbors of v (the degree of vertex v). Then the clustering coefficient of 
vertex f , c^,, is given by: 

26u 

dy{d^ - 1) 

Note that total number of possible connections among neighbors of v is 

2 

Therefore, the clustering coefficient is essentially the probability that two 
neighbors of v are connected. For a graph G with n vertices, the average 
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clustering coefficient, c, is defined as: 

1 "■ 



c = - y c„ 
n 

v=l 



3.6. Generating a Related Random Graph for Small-World Analysis. 

In order to analyze our networks for small-world properties, it was useful 
to compare these graphs to those of similar networks with randomly as- 
signed edges. Small- world networks are nearly as well-connected as random 
graphs, but possess a surprisingly well-localized structure. We developed 
the following algorithm for this process: 

For a graph G with n vertices, let k be the number of edges on G. There- 
fore the average number of edges per vertex is k/n. Furthermore, the proba- 
bility that any two random vertices are connected, p, is given by the number 
of existing connections divided by the total possible connections: 

k 2k 

P 



"("-^ n(n - 11 
2 ^ ^ 



Next we generate a random graph, Rand(G), with n vertices and a p proba- 



bility of connection between two vertices (See Methods 3.1 1. We then com- 
pute c and / for G and Rand{G). 



G has small- world network properties [33] if: 

1. 1{G) > l{Rand{G)) 
and 

2. c{G) > > c{Rand{G)) where c{Rand{G)) k. p 



3.7. The Neighborhood of a Graph. On a graph G with n vertices, the 
neighborhood, H{m) where m is a positive integer, is useful in analyzing 
small- world networks. A neighborhood of size m around vertex v is the set 
of all vertices that can be reached from v mm steps or less. We shall refer 
to the number of vertices reachable in m steps or less as Hy(m). To get a 
sense of the global neighborhood size on the graph, we can average H^im) 
over each vertex v in G. This gives us a global H{m) for a given m: 

1 " 

Him) = -y^HJm) 
n ^-^ 

v=l 

It is clear that H{m) is non-decreasing and as long as G is complete, H{m) 
achieves its maximum, n, for finite m. 
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3.8. Graph Rewiring. The rewiring principle can be rigorously explained 
as follows. We first number each vertex in G from 1 to |G| (where n— \G\,n 
being the total number of vertices). If there is a connection between vertices 
u and V in G, we generate a random number between and 1 . If this random 
number is less than a given probability p, then the connection will be rewired. 
Without loss of generality assume u < v. We then fix the connection to 
vertex u, and move the connection to another vertex, k, such that u and k are 
now connected whereas they were not previously. 

3.9. Graph Energy. For a graph G = {V, E) where V is the set of vertices 
and E is the set of edges, one can define an arbitrary scalar function u : 
V — > M. Consequently, one can then define the energy of u associated with 
the graph. The energy of u associated with G, £{u) is defined as: 

£{u)= J2{u{x)-u{y)r. 

x,yeE 

We analyzed the energies of the Laplacian matrix eigenfunctions, thus u — 

3.10. Spacial Variance. In order to discuss spacial variance, we must first 
define the resistance between two vertices on a graph. Let G — (V, E) be a 
graph and x,y eV . Then the resistance between x and y, d{x, y), is given 
by: 

d{x,y) = £{h{x,y)) ^ 

Where h{x, y) is a harmonic function defined as follows: 

Let G — (y, E) be a graph and x,y & V. Then the harmonic function 
corresponding to {x, y) is a scalar function h{x, y) : F — > R such that: 

1. h{x, y){x) = 

2. h{x,y){y) = 1 

3. 8{hY^, where h is an arbitrary scalar function on V, is maximized at h{x, y). 

Finding the harmonic function is equivalent to finding a vector h such that 
Lh = z, where z is a vector whose entries are all except for those entries 
corresponding to x and y. This is analogous to what "harmonic" means in 
Euclidean space. This changes the maximization problem in condition 3. to 
solving a system of linear equations. 

Using these we can now define the spacial variance of a graph. Again, let 
G — {V, E) be a graph with n vertices and be a scalar function on V. Let 
7 be a constant. Then the 7*'^ spacial variance of u over G is given by: 

Var^{u) = - d{x,yy{u{x) - u{y)f 

x,yeE 
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In this paper, the spacial variances of eigenfunctions of Laplacian matrices 
were evaluated at 7 = 1 and then analyzed. 
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